Spin Order due to Orbital Fluctuations: Cubic Vanadates 
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We investigate the highly frustrated spin and orbital superexchange interactions in cubic vana- 
dates. The fluctuations of t2g orbitals trigger a novel mechanism of ferromagnetic interactions 
between spins S = 1 of V"^"*" ions along one of the cubic directions which operates already in the 
absence of Hund's rule exchange Jh, and leads to the C-type antiferromagnetic phase in LaVOs. 
The Jahn- Teller effect can stabilize the orbital ordering and the G-type antiferromagnetic phase at 
low temperatures, but large entropy due to orbital fluctuations favors again the C-phase at higher 
temperatures, as observed in YVO3. 
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Large Coulomb interactions play a crucial role in tran- 
sition metal oxides, and are responsible for the collective 
behavior of strongly correlated d electrons which localize 
in Mott-Hubbard (or charge-transfer) insulators . Such 
localized electrons may occupy degenerate orbital states 
which makes it necessary to consider orbital degrees of 
freedom at equal footing with electron spins, and leads to 
the effective (superexchange) spin-orbital models to de- 
scribe the low-energy physics . A remarkable feature 
of these models is that the superexchange interaction is 
highly frustrated on a cubic lattice, which was recognized 
as the origin of novel quantum effects in transition metal 
oxides In case of Cg orbital systems this frustration 
is likely removed by orbital order due to order-out-of- 
disorder mechanism, which maximizes the energy gain 
from quantum spin fluctuations 1^. Moreover, quantum 
effects among Cg orbitals are largely suppressed by the 
Jahn- Teller (JT) effect in real systems, which together 
with superexchange often leads to structural phase tran- 
sitions accompanied by a certain ordering of occupied or- 
bitals, supporting particular magnetic structures. Some 
well known examples are systems with degenerate eg or- 
bitals filled either by one hole (KCuFs), or by one elec- 
tron (LaMnOs), which order antiferromagnetically well 
below the structural transition. 

The transition metal oxides with partly filled t2g or- 
bitals exhibit different and more interesting phenomena. 
This occurs due to the relative weakness of the JT cou- 
pling in this case, and due to the higher degeneracy and 
additional symmetry of t2g orbitals [Q. As a result, the 
orbitals may form the coherent orbital-liquid ground state 
stabilized by quantum effects, as observed in the spin 
S = 1/2 Mott-insulator LaTiOs It is puzzling what 
happens when the t2g orbitals are filled by two electrons, 
as in vanadium oxides. On one hand, the occupied t2g 
orbitals are known to order in non-cubic vanadium com- 
pounds, such as LiV02 and V2O3 In fact, the 
first spin-orbital model for V2O3 with spins s — 1/2 was 
proposed over twenty years ago fO], but later it was re- 



alized that Jh at Y^^{d^) ions is large j|T2|, and the rele- 
vant model has to involve 5 = 1 spins ]10| ]. On the other 
hand, the situation in cubic systems might be very differ- 
ent as all the bonds are a priori magnetically equivalent, 
and quantum fluctuations among orbitals are expected 
to play an important role in this case. 

In this Letter we derive and investigate the spin-orbital 
model for cubic vanadates: LaVOa and YVO3. The 
magnetic order in LaVOs is C-type [ferromagnetic chains 
along c-axis which stagger withi n (a, b) planes], with the 
Neel temperature Tjv = 140 K ||l|,|lj, while it is stag- 
gered in all three directions (G-type) in YVO3 at T < 77 
K and C-type at higher temperatures 77 < T < 114 K 
Jl^-p^. The C-phase is particularly surprising as aris- 
ing from a practically undistorted structure of LaVOa at 
T > Tn jl^. Recent Hartree-Fock studies have shown 
that indeed C- and G-phase are energetically close [ pT| . 
In order to understand the microscopic origin of their 
competition we consider the regime of large ?7, and ad- 
dress below the following questions: (i) can the superex- 
change interactions alone explain why the ferromagnetic 
(FM) and antiferromagnetic (AF) interactions coexist in 
LaVOa in spite of a practically ideal cubic structure at 
T > Tpf with almost equal V-V bonds; (ii) why the struc- 
tural transition in LaVOa occurs only below the magnetic 
transition; and (iii) why is the G-type AF order stable 
in the low-temperature phase of YVO3, while the C-type 
order wins at higher temperatures? 

We start with a Mott-insulator picture of cubic vana- 
dites, consistent with the large value of on-site intraor- 
bital Coulomb element U ~ 4.5 eV [Q, and with the 
results of electronic-structure calculations [fsl. Due to 
the Hund's coupling Jh — 0.68 eV @ the ions are 
in triplet configuration •^T2. Each t2g orbital is orthogo- 
nal to one cubic axis. For instance, yz is orthogonal to 
a axis and will be labelled as a, while zx and xy will be 
labelled as b and c, respectively. The electron densities 
at V'^"'" ions satisfy a local constraint, rija +nn, + nic = 2. 

The superexchange interactions between 5 = 1 spins 
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arise from the virtual excitations dfd3 dfd] on a given 
bond (ij), with the hopping t allowed only between two 
out of three t2g orbitals. The excited state may be ei- 
ther a high-spin state, or one of three low-spin states: 
"^E, ^Ti or Jll. When the second order processes 
d'^d'j dfd] ^d^d'j are analyzed, one has to project 
the d^ idf) configuration generated after an individual 
hopping process on the above eigenstates (^T2 ground 
state). This leads to the spin-orbital Hamiltonian, 

^ = ^E E [{s^-s, + i)j[]^ + k\]^]^ (1) 

T (y'>ll7 

where the orbital operators j'l^'^ and k\'^'^ follow from the 
processes active on the bond (ij) |1 7, where j = a,b, c: 

jI]^ = \[{l + 2r/i?) (t, • f, + \n,n^) 

~ W [t^t] + \n,n^) - i,7i?(n, + n,)] ^''^ , (2) 
k'i]^ = ViR {^i ■ fj + jUifij) + Tjr {t^t^ + in^Tij-) 

- i(l + 77^)(n,+7i,)]^''\ (3) 

and J = U'^/U. The coefficients i? = 1/(1 - 877) and 
r = 1/(1 -|- 2?]) originate from the multiplet structure of 
the t\g excited states via 77 = Jh/U |i^. The operators 
Ti = {t^ ,Tf ,tI} are defined in the orbital pseudospin 
subspace spanned by two orbital flavours which are ac- 
tive along a given direction 7. For instance, for a bond 
(ij) II c, the interactions follow from the electron hop- 
ping between the pairs of a and b orbitals, and may be 
expressed by Schwinger bosons: = a|6j, t~ = bla^, 

= 5("« - "if,)> and ' = n„ + Ti^ft, and therefore: 

2 {n ■ Tj+ jniUj^"^ = {n.anja + af^^^^aj) + (a <^ b), 
2 (t/ t/ -f iniUj) = n^anja + ^J»b"j6- (4) 
Consider first the interactions in the Jh —>■ limit: 

where a constant energy of —2 J per ion is neglected. 
It is straighforward to understand why the interactions 
in this limit turn out to have the same structure as in 
LaTiOs 0, where for s = 1/2 spins of Ti^+ ions one finds 
instead {si ■ Sj + j). In fact, the spin interactions follow 
entirely from the Pauli principle, as the terms oc Si ■ Sj 
due to the high-spin and low-spin '^E states, which 
involve d'^{abc} configurations, cancel each other. 

A remarkable feature of the t2g superexchange in Eq. 
is that every bond is represented by two equivalent 
orbitals giving a SU(2) symmetric structure (ji ■ fj + 
jn^nj)^'^') of the orbital part. Depending on the type 
of orbital correlations this may result in a spin coupling 



constant of either sign. This important property resem- 
bles that of one-dimensional (ID) SU(4) model |0. The 
present problem is however more involved since there are 
three t2g flavours in a cubic crystal, and SU(2) orbital 
correlations among two of them along a particular direc- 
tion will necessarily frustrate those correlations in other 
directions. One may also notice a certain analogy with 
the models of valence bond solids |^ : Actually, a large 
orbital moment L = 1 of t2g states is formally decom- 
posed in Eq. (^) onto pseudospins one-half, active on 
different bonds. The analogy is again only partial since 
independent rotations within orbital doublets active on 
different bonds are not allowed here by construction, and 
thus the formation of orbital singlets in all directions si- 
multaneously is impossible. 

Another key observation is the difference between the 
interactions derived for the pairs of Ti^+ and V'^"'" ions: 
In the d^ configuration spin s = 1/2 is small, and the idea 
of composite spin-orbital resonance, in analogy to SU(4) 
excitations |^^, helped to resolve the orbital frustration 
problem 0. In that case the superexchange is best op- 
timized by the spin-orbital resonance in all three direc- 
tions, and the orbitals form a three-dimensional quan- 
tum liquid which coexists with weak spin order of G-type. 
This mechanism is however suppressed in the present case 
of large spin S* = 1 at ions, and the quantum energy 
can be gained mainly from the orbital sector. This sug- 
gests that a particular classical spin configuration may be 
picked up which maximizes the energy gain from orbital 
fluctuations. Indeed, orbital singlets (with nia+nn, = 1) 
may form on the bonds parallel to c axis, exploiting fully 
the SU(2) symmetry of the orbital interactions in one 
direction, while the second electron occupies the third 
t2g orbital (n^c = 1), controlling spin interactions in the 
(a, b) planes. 

In order to understand why orbital fluctuations sup- 
port the C-AF type spin order, it is instructive to start 
with a single bond along c-axis. A crucial observation 
is that the lowest energy of — J/2 is obtained when the 
spins are ferromagnetic, and the orbitals a and b form a 
singlet, with {fi-Tj)'^'^^ = — | |2^. Thus, one finds a novel 
mechanism of ferromagnetic interactions which operates 
due to local fluctuations of a and b orbitals. At the same 
time, the orbital resonance on the bonds in (a, h) planes 
is blocked, = 1 , and the superexchange is AF 

due to the excitations to ^Ti and '^T2 states with doubly 
occupied c orbitals. Such an electron distribution and the 
formation of quasi-lD orbital pseudospin chains supports 
FM spin order along c axis in C-phase, and could be sta- 
ble only at low temperature when a coherent spin state 
breaking the cubic symmetry is formed as well. The on- 
set of the magnetic order which coexists with such orbital 
fluctuations explains also why a structural transition is 
here induced by this coherent electronic state. 

We compared the ground state energies of C- and G- 
phase using the spin- wave theory for the spin part, while 
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either the exact Bethe ansatz result, or the Gaussian fluc- 
tuations around the ordered state, were used for the or- 
bital part. The exchange constants within (a, b) planes 
(Jab > 0) and along c axis (Jc < 0) determine the spin 
waves. They follow from Eqs. (in units of J): 



Jc 



1 - ri{R+r) + (1 + 2t]R - rir){nianja) 



Ac) 



J _ 1 

•Jab A 



\{b) 



(6) 



and are given by orbital correlations. Their values at 
= are easily obtained from the Bethe ansatz result for 
a ID Heisenberg chain of disordered a and h orbitals (in 
this case {nianja)^''^ = |): Jab = 0.313 and Jc = —0.097. 
Although the FM interaction along c-axis Jc is weaker in 
this limit, it gives a considerable energy gain of ~ 0.2 J 
over the G-AF order. It is further enhanced at Jh > 
by a mechanism similar to that known from cuprates 
and manganites as the high-spin '*^2 state lies by 

3Jh — 2.0 eV below its low-spin "^E counterpart. This 
splitting modifies the ID orbital-wave (OW) spectrum. 



= VA2+i?2(l-cos2 k), 



(7) 



and the gap A = {ri{R + r)[2R + r]{R + r)]}i/2 

opens. 

Using the spin-wave theory we determined the orbital 
correlations in Eqs. (||). As a result, one finds increas- 
ing (decreasing) FM (AF) exchange constants with in- 
creasing Jh (Fig. |l|). Taking a representative value of 
the hopping integral t — 0.2 eV gives J = 35.6 meV 
which leads to the exchange constants in the C-phase ob- 
tained for a realistic ratio Jh/U = 0.15: Jab — 7.1 and 
Jc — —9.3 meV. These values are in the expected range, 
taking the Neel temperature Tjv — 140 K of the C-phase 
in LaVOa. We emphasize that the orbital quantum fluc- 
tuations play here a dominating role and the well known 
Hund's mechanism due to Jh alone would not suffice to 
obtain \Jc\ > Jab, giving instead Jc ~ —4.4 meV. 
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FIG. 1. Exchange interactions Jc and Jab (j^) in C-phase 
(in units of J) as functions of = Jh/U. The inset indicates 
the local fluctuation between a (= yz) and b {= zx) orbitals 
due to singlet formation along the c axis. 

Next we consider the reasons for the stability of the 
G-phase in YVO3. Unlike LaVOa with almost equal V- 
V bonds [Q, this compound crystafizes in the distorted 



structure [^,^, indicating that the JT effect plays a 
significant role. It was suggested that energy may be 
gained due to C-type orbital ordering, with a and b or- 
bitals staggered in (a, b) planes and repeated along c-axis, 
while riic — 1 pq-p^ . Such ordering can be promoted by 
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and competes with the orbital disorder. This behavior is 
remarkably different from the eg systems, where the JT 
effect and superexchange support each other, inducing 
orbital ordering While U > causes orbital split- 

ting by 4U and thus lowers the energy of the G-phase 
(Eq), it has little effect on the energy of the C-phase 
(Ec). The energy difference is given by 



AE^Ec-EgO:^V- \ri{m 



5E, 



OW, 



(9) 



where SEqw > is the energy contribution due to quan- 
tum fluctuations of t2g orbitals. Large spins S = 1 are 
almost classical and their fluctuations could be neglected. 
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FIG. 2. Free energies J-{T) (in units of J) of: G-phase ob- 
tained with V = 0.65 J (full line), and C-phase for 77 = 0.05, 
0.10 and 0.15 (dashed lines) , as functions of temperature T/J. 

Orbital excitations are quite different in both AF 
phases: while the gap A is small and grows cx in C- 
phase (0), a large gap ~ 4U opens in the OW spectrum 
of G-phase, uj^ — AV + rjRcosk. Thus, both the larger 
quantum fluctuations and additional (classical) energy 
gain due to finite Jh in the C-phase have to be overbal- 
anced by the JT energy cx U in order to stabilize the 
G-AF order at T ^ 0. However, the G-phase may be 
destabilized at finite T by larger orbital entropy of the 
C-phase. Indeed, taking V — 0.65J and 77 = 0.15, the 
free energy, = {H) — TS, with the entropy S deter- 
mined by orbital excitations, gives a transition from G- 
to C-phase around T* ~ 0.8 J (Fig. |). While this be- 
havior reproduces qualitatively the first order transition 
observed in YVO3 [ p^p^ , its quantitative description re- 
quires a careful consideration of lattice and spin entropy 
contributions to the free energy T. These effects are ex- 
pected to reduce the transition temperature T* down to 
experimental values. 
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The exchange constants in the G-phase are anisotropic 
which can be understood by analyzing the superexchange 
expressions (^). The transitions to the low-spin ^Ti and 
states occur in all three directions between pairs of 
occupied orbitals of the same kind, and give the leading 
AF contribution cx (1 — rjr). The excitations of d'^{abc} 
configurations occur in addition on the bonds in the (a, b) 
planes, and reduce the AF coupling Jab by rjR, giving 
Jab < Jc- Including in addition a relativistic spin-orbit 
coupling oc A (~ 18 meV [E2|), we find for the G-phase: 
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Je = i(l-77r)-i(l + 277i?-r7r)A2, 
Jab = Ul-vR- Vr) + i(l + 277i? - 



(10) 



where A = A/AV. The spin-orbit coupling enhances (re- 
duces) the effective superexchange in (a, b) planes (along 
c-axis), as shown in Fig. ^. For example, taking a value 
of A/4V = 0.2 which is believed to be close to realistic for 
YVO3 , we find Jab — 5.9 and Jc — 6.9 meV. These values 
lie in the expected range for the G-phase of YVO3. This 
also demonstrates that the magnetic structure and the 
spin-wave spectrum are completely different depending 
on the orbital state — the exchange constant Jc which is 
FM in LaVOa may become even the strongest AF bond 
when the orbitals have ordered, as in YVO3, and the JT 
splitting (H) dominates over the spin-orbit coupling. 
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FIG. 3. Exchange interactions Jc and Jat ( |ic| ) in G-phase 
(in units of J) as functions of spin-orbit coupling A/AV, for 
•q = 0.15. The inset shows a and b orbitals, which are stag- 
gered in (a, b) planes, and repeat themselves along c axis. 

Summarizing, strong t2g orbital fluctuations in a half- 
filled system of yz and zx orbitals in cubic vanadites 
lead to a new mechanism of ferromagnetic superexchange 
which stabiHzes the C-phase in first undistorted LaVOs, 
and the structural transition follows. The JT effect op- 
poses the superexchange and can stabilize the G-phase 
with orbital ordering but only at low temperatures, as 
the fluctuations of t2g orbitals release high entropy, and 
are thus responsible for the transition from the orbital 
ordered G-phase to the orbital disordered C-phase, ob- 
served in YVO3. 
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